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1. INTRODUCTION 


The serious study of the phenomenon commonly known as the 
pressure broadening of spectral lines is now over three-quarters of a 
century old. It was in 1895 that Michelson made the first important 
contribution to this subject, and identified the mechanism of pressure 
broadening as the sudden interruption of the vibrations of the atomic 
radiation by collisions with neighbouring atoms. A simple Fourier 
analysis then shows that the line width is of the order 1/7 where 7 is the 
mean time between collisions. This fundamental idea was followed up 
by Lorentz (1906) who concentrated his attention on absorption rather 
than emission and showed that the effect of random sudden collisions 
is equivalent to an additional damping force, proportional to velocity, 
exerted on the atomic oscillators. The line profile then has the well- 
known Lorentzian shape. 

A quite different approach to pressure broadening was initiated by 
Debye (1919) and Holtsmark (1919) in the context of perturbation of 
the radiating atom by charged particles. Here the environment of the 
radiator is considered to exert a slowly changing influence on the 
frequency emitted, and the intensity profile of the spectral line is given 
by the probability of occurrence of the configurations of perturbers 
which give the various frequency displacements. The application of this 
idea to broadening by neutral particles was made by Jablonski (1931) 
and, in more detail, by Margenau (1932). This type of treatment of the 
problem is known as quasi-static (or statistical) theory. 

Important progress was made by Weisskopf (1932, 1933), who 
perceived that the nature of the interruption of the atomic oscillations 
visualized by Michelson and Lorentz was essentially an interruption of 
phase. If the phase change induced by the collision is greater than a 
certain critical value, which Weisskopf took to be | radian, then the 
wave train is effectively terminated, and it is the mean time between 

143 
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such “strong” collisions which determines the width of the resulting 
Lorentzian profile. If the interaction potential between the perturbing 
and radiating atoms is known, it is possible to calculate the line width. 
Weisskopf further showed that a simple quantum mechanical calcula- 
tion based on the same physical idea leads, in the WKB approximation, 
to the same result. 

In another section of his wide-ranging review, Weisskopf (1933) 
discussed the relationship between the two apparently different 
approaches to pressure broadening. The interruption, or impact, 
approach entirely neglects the radiation emitted during a collision, 
while the quasi-static approach takes into account radiation emitted at 
all times but does not make a Fourier analysis and consequently is a low 
velocity limit of a more general theory. 

Since the work of Weisskopf, progress in the theoretical understand- 
ing of pressure broadening caused by neutral atoms has developed along 
four main lines: 


(1) The impact theory for classical atomic oscillators has been com- 
pletely developed, taking account of collisions which induce small phase 
changes (“‘weak’’ collisions) as well as those which induce large phase 
changes (“‘strong”’ collisions). The former collisions give rise to a shift 
of the spectral line, which, however, remains Lorentzian in shape (Lenz, 
1933; Lindholm, 1945; and Foley, 1946). 

(2) Attempts have been made to develop a quite general classical 
theory which includes the impact and quasi-static approximations as 
limiting cases. Lindholm (1945) achieved an approximate solution to 
this problem. Other attempts are referred to in the earlier review articles 
mentioned below, and recent attempts are discussed in section 2C of 
this article. 

(3) The quantum mechanical theory of impact broadening has been 
the subject of many contributions in recent years. An early attempt to 
establish the general formalism by Spitzer (1940) was followed by 
further elucidation by Anderson (1949), Baranger (1962) and Griem 
(1964), and from a different point of view by Fano (1961). Resonance 
broadening (where the perturbing and radiating atoms are of the same 
species and can exchange excitation energy in a collision) has attracted 
a great deal of attention, and important contributions have been made 
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by Byron and Foley (1964), Omont (1966), Mead (1960) and Ali and 
Griem (1965, 1966). Some of these matters form the subject of section 
2B of this review. 

(4) In the field of “foreign” gas broadening, with which this article 
is mainly concerned, developments have largely centred on the question 
of the interatomic potential. The use of the Lennard-Jones potential 
has been explored in the impact approximation by Behmenburg (1964) 
and by Hindmarsh and co-workers (1963, 1967), and in the quasi- 
static approximation by Bergeon, Robin and Vodar (1952) and Bergeon 
(1954) and by Farr and Hindmarsh (1968). Section 2C of the present 
work is devoted to these developments. 

The most important of the previous reviews of this subject is that of 
Ch’en and Takeo (1957). This provides a comprehensive account of 
both experimental and theoretical work up to that date. In the present 
article we shall mention work earlier than 1957 only where it is necessary 
for understanding the later developments. Other useful reviews and 
books are by Sobelman (1954), Unsdld (1955), Traving (1960) and 
Breene (1957, 1961). The present work does not aim at total compre- 
hensiveness, but rather at establishing the present position achieved by 
theoretical and experimental investigation. 


2. THEORY 
A. Classical Theory of Collision Broadening 
1. General formulation 


The classical model represents an atomic radiator as an elastically 
bound oscillating electron, and the problem of collision broadening is 
to calculate the spectrum of such an oscillator when it is disturbed by 
collisions with the other atoms in its environment. We assume that the 
collisions do not change the amplitude of the oscillator, but only the 
phase. This is equivalent to assuming that the transition probability is 
unchanged by the collision. This seems on the whole to be confirmed by 
observation, at least at low pressures, although it has been suggested 
(e.g. Jefimenko, 1964) that in some cases it is incorrect. The evidence 
is not altogether clear, but all existing theories involve this assumption. 
The treatment of the classical problem which follows is largely based 
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on that of van Trigt (1966) who even at that late date found that 
several points required clarification. The displacement of the oscillator 
as a function of time is given by 


t 


f(@Z) ea.Re jexp {i i w(t’) art t>0O (1) 


where w(t’) is the instantaneous angular frequency of the oscillation 
at time ¢’ and Re stands for “‘real part of”. The undisturbed frequency 
will be represented by wo, and the oscillation begins at t = 0. The effect 
of natural damping is ignored. It is, in fact, often negligible compared 
with collisional effects, but in the impact approximation the effects of 
natural and collisional damping are additive, as Lorentz (1906) showed. 
The intensity profile J (w) of the light emitted by the oscillator is then 
given by the squared modulus of the complex Fourier transform of f(t): 


I(w) = | { {re | exp \! few art} exp {iwt} dt *. (2) 
0 } 


A more manageable, and entirely equivalent, general expression is 
found by first writing the autocorrelation function of the disturbance 
f(t). Putting 


t 


| w(t) dt! = wot + ft) 


10) 
where 7(t) is the departure in phase of the perturbed oscillation from 
the unperturbed oscillation, the autocorrelation function G(s) is: 
@(s) _ [f(t) f(t = Sy —_ 4[Re fertmot ei os eine nay 
+ Re {ei@os ete ray 


where [ ], indicates the average value over all values of t. Now the 
first term in the expression for ®(s) is zero provided that no significant 
disturbance in phase occurs, on the average, within one cycle of oscilla- 
tion. For optical spectra, where the linewidth is a small fraction of the 
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undisturbed frequency, this is a very good approximation, but it 
should be noticed that the results obtained with its use cannot be 
extrapolated to give the correct zero frequency results. For our purposes, 
then, we may write 


(5) = H{Re [olor olM-9]}, 
See er re ee 
=(6) ahi *hs) (3) 
where n(t,s) = y(t + s) — 7(Z). 


Finally we recall that the intensity profile is the complex Fourier trans- 
form of the autocorrelation function of the disturbance: 


I(w) = | Ws) eS ds + | P*(syerds: (4) 


No further progress can be made without evaluating the time averages, 
and to accomplish this approximations must be made. 


2. The impact approximation 
We first write down an expression for A[¥(s)e~'®°5], a change, not 
necessarily small, in ¥(s)e~'®°5 corresponding to a change As in s: 


A[Y(s)e“'**"] = dle"? (e'4" — 1)], 
= dei"), [e'4” — 1], 
= P(e" [ot — 1], 
where An = 7(t,s + As) — r(t,9) 


and the second line follows because the changes Ay occurring in an 
oscillator perturbed by random collisions are statistically independent 
of 7, the accumulated phase disturbance at time s. We next suppose 
that the phase changes Ay are caused by successive collisions. This is 
known as the impact approximation and its range of validity will be 
discussed later in this section. The number of collisions with impact 
parameter p which occur in the time interval As is Nid 27zpdp|As| where 
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N is the number density of perturbing particles and the further assump- 


tion has been made that the perturbing and radiating atoms move in 
straight lines with mean relative velocity 6. We then have 


[e'4" — 1], = 27Nd| As | | [e'"™ — 1] pdp 
oO 


where 7(p) is the phase change caused by a collision of impact para- 
meter p. We write the integral in the form 


[e'4" — 1], = — Nod| As | (c, — io;) 
where o == Jr | [1 — cos 7{p)] pdp (5) 
0 
and o; = 27 | sin 7 (p) pdp. (6) 


c0) 
It follows that 
A[P(s)e7'2] 
——————— = — No| A AE 
sje ioe Nd| As | (o, — io;). 


Since the average over a large number of collisions is iavolved, this 
last expression can be treated as a differential equation. Integration then 
gives 


Y(s) = constant X exp {[iw) — Nd (c, — io,)]| 5 |}. (7) 


Thus the spectrum is, from eqns. (4) and (7), 


es) SON Ge. —————— 
(@> + w + Noo)? + (Nio,)? 


+ Nord ein Jaieiatea sla ; 8 
(wo — w + Noo)? + anit (8) 
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In optical spectra the first of the two terms in the final expression for 
I (w) is negligible compared with the second. Neglecting the first term, 
and normalizing so that 


fo 0) 


| I(w) dw = 1, 
0 

we obtain 

Noo,/7 


(| 
3) (wo — w + Néo;)? + (Nio,)” 


(9) 
The spectral line therefore has a Lorentzian shape, but the whole line 
is shifted with respect to the undisturbed line centre at wo. It is clear 
that the quantity o, — io; can be looked on as a complex cross-section, 
and that the broadening of the line is associated with the real part o, 
of this cross-section and the shift with the imaginary part —o;. It is 
also easy to see from eqns. (5) and (6) that, since 7(p) increases rapidly 
as p decreases, the contribution to o; arises predominantly from small 
values of 7 while that to o, arises predominantly from large values of 7. 
This corresponds to the original physical idea of Weisskopf (1933) that 
the broadening arises essentially for collisions which cause a phase 
change greater than about | radian. 

All that remains to press this classical calculation to a numerical 
conclusion is to evaluate 7(p) and then perform the integrations in 
eqns. (5) and (6). If the interatomic potentials of the upper and lower 
states of the radiating atom and the perturbing atom (taken to be in 
the ground state) are respectively V,(R) and V,(R), where R is the 
distance between the two atoms, then 


+0 


| t. 1 
no) =; | WAR) — rarer => [CR at 


—0 


where f represents the time measured from the time of closest approach 
of the colliding atoms. If we assume that the atoms move in straight 
paths with velocity 6, then, 


R2 =< p? + D7 2 
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so that 


mete: ” Y(R)RAR 
HO Fy (R? — pave (10) 


The question of a suitable form for the functions V(R) will be taken up 
in section 2C. 

We now turn to an examination of the validity of the important 
approximations which have been made. 

First we consider the impact approximation itself. Although much 
has been written about this (see the reviews mentioned in the Introduc- 
tion), the essential idea is very simple. The collisions are supposed to 
occur successively and to occupy a time small compared with the time 
between collisions. The radiation emitted during the collisions may 
then be neglected. Now the duration of a collision is r, ~ p/i and the 
mean time between collisions which cause broadening is 7 = (Nio,)7}. 
We compare the intensity per unit frequency interval from 7, with that 
from 7 at a given frequency displacement (w) — w) from the line centre. 
Remembering that the former contribution is spread over an interval 
~7, *, the ratio of the two contributions is approximately 

(p[i)? (wy — w)? _ p(y — w)? 
(Nic,)~1 Noo, 5? ; 
If the impact approximation is to be valid this ratio must be < 1. 
Hence 


(wo — w)? < 0?/mp?. (11) 


But the value of w) — w associated with a certain value of p is deter- 
mined by the interatomic potential. We put 


wo —w = — V(p)/h 
and write the solution of this equation as 
p = Ww — w). 


The inequality which defines the region of validity of the impact approxi- 
mation may then be written 


(wo — w)?W? < b?/z. 
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If |V(p)| is a single-valued function of p which decreases more rapidly 
than linearly with p, then (w) — w)?W is an increasing function of 
|w) — w|. The impact approximation is therefore seen to be valid at 
small frequency displacements and to be essentially a high velocity 
approximation. If V(p) is not a single-valued function, and in practice 
it never is (but has the characteristic shape of an interatomic potential), 
the impact region is more difficult to define. But in most cases, 
the general conclusion already mentioned is valid. If we take the 


van der Waals interaction, V(p) = — C,/p®, giving W(wo — w) = 
(C,/h)*/°/(wo — w)*/°, the inequality (5) becomes 
D 6/5 ( h Jes 
jo — | < (2) C. é (12) 


For neutral atom pairs Cs ~ 107°’ erg cm®, and for many experiments 
6 ~ 10° cm sec~*. With these values, we find the impact region defined 
as 

fier —tn | <3) 10** rad sect, ==-5.cm=1, 


Although values of C, and é can, of course, vary over a certain range, 
the impact region is mostly defined by a boundary of a few wavenumbers. 
It is of some consequence to notice that this boundary occurs on one 
wing only of the line. If Cg is positive (the usual case), then wo — w is 
positive and the boundary refers to the low frequency wing of the line. 
If the potential is entirely described by the van der Waals interaction, 
the impact approximation is, according to the above account, valid for 
all positions on the high frequency wing of the line: all radiation emitted 
during the collisions is displaced to longer wavelengths. In fact, for 
very close encounters the “repulsive” part of the potential curve causes 
radiation to be emitted on the high frequency side of the line centre, 
so that in practice an impact boundary arises on that side also. Taking 
the Lennard-Jones form V(p) = C,,/p'* — C,/p® it turns out, with 
typical values for the constants, that the impact boundary on the low 
frequency side is not usually affected by the presence of the R~'? term, 
and that on the high frequency side the boundary is determined only by 
the repulsive term. With C,, = 10~1°* erg cm?” and 6 = 10° cm sec™?, 
we find 


agi w) ae 10 vad sec" * = T5'em=** 
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Secondly we examine the validity of the assumption that the atoms 
move in straight lines. A qualitative calculation can be made to estimate 
the range of applicability of the straight path approximation. 

The orbit will be significantly curved if the impulse of the force F 
exerted between the colliding atoms becomes comparable with the 
momentum of either of them. The relevant ratio is therefore 


_ frat @V/AR)a=, 6 bee = Zs 
~ Mob Mo MB? 6 


Li mid 
D ake p 
where M is the mass of the smaller atom involved. Now the phase 
change, calculated for straight paths, for a Lennard-Jones potential is 
(Hindmarsh, Petford and Smith, 1967) 


7 Pact eh er (13) 


If, at the critical value p, of p, 7 is determined essentially by the second 
term, then p is also so determined. We then have, since 7 = 1 at p,, 
De Os C,/ho 
and 
6h 6h°/> 
pr = ee 
Méip. . Mec 
on substituting typical values. Path curvature is therefore small. 
If Cg is small, p, is essentially determined by the R~/? term. In this 
case we find p ~ 1 with C,, = 107112 J cm'?. We therefore expect 
path curvature to have a marked effect. A fundamental point which 
arises here is how to combine the potentials for the two levels of the 
radiator which are involved in the transition, since the classical path 
will not be the same for both levels. This has been discussed by Fiutak 
and Czuchaj (1970) who point out that the condition p < 1 for straight 
paths is essentially the same as the condition for the pertuber motion 
to be described classically at all (see section 2B 1). Although Schuller 
and Vodar (1960) have calculated the effect of curvature of the orbits 
for a Sutherland type of potential, and find significant effects in certain 
circumstances, a greater effort is probably not worth while unless a 
full quantum mechanical treatment is used. 
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Thirdly we mention the problem of velocity averaging. In most of 
the literature velocity averaging is achieved by the simple device of 
writing 3, the mean thermal relative speed of the radiating and perturb- 
ing atoms wherever the velocity enters the calculation. This is clearly 
not correct. It is a relatively simple matter to average the phase shift 
expressions over a Maxwellian velocity distribution, and this in fact 
makes only a small difference to the numerical results. The real diffi- 
culty is that if rigorously correct velocity averaging is to be carried out, 
the thermal (Doppler) broadening of the spectral line cannot be treated 
separately from the collision broadening as is usually done. This 
question has been studied by Gerston and Foley (1968), who refer 
to the earlier work of Dicke (1953) and Wittke and Dicke (1956) on 
““pressure narrowing”. The essential point is that if the mean free path 
of the radiating atoms for collisions which change their direction be- 
comes much smaller than the wavelength of the radiation, the Doppler 
shift induced in one portion of a cycle of the radiation tends to cancel 
the shift from another portion. The effective Doppler width is then 
much reduced, and provided the ordinary collision broadening is not 
too great, a net narrowing with increasing pressure can occur over a 
certain pressure range. This is not usually the case in optical spectra, 
but at much lower frequencies the phenomenon is well known. 


3. The quasi-static approximation 

The quasi-static approximation takes into account radiation emitted 
at all times, but does not make a Fourier analysis. The intensity at an 
angular frequency w is thus proportional to the fraction of the total time 
which the atom spends in radiating at that frequency. If the instantan- 
eous frequency is represented by w = wy + Aw then the quasi-static 
theory gives 


I(Aw) oc P(Aw) 
where P(Aw) is the probability of occurrence of the frequency displace- 
ment Aw. We now see how this follows from the general theory already 


outlined. We begin once again with the autocorrelation function (s) 
of the disturbance and work as before with the function (s) defined by 


O(s) = Y(s) + ¥*(s). 
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Then 


Ys) =— feiees (aiden Be 


The essence of the quasi-static approximation is to assume that 7(t,s) 
is a slowly varying function of s, so that in the time that 7(¢,s) has 
changed significantly, e’@o° has undergone a large number of oscilla- 
tions. We then expand 7(t,s) in a Taylor series about f: 


d 
n(t,s) = ft +3) — a= s+... =sho +. 


since Aw = dy/dt. 
Neglecting all but the first term in the expansion, we have 


Ys) = ef@osfelsAo] 
If, then, P(Aw) is the probability of occurrence of Aw, we may write 


+ 00 
Ys) = | els gi8A@ P(Aw) dAw 


—-o 


and the line profile is 


+0 +0 
I(w) = | | ei(@+ao)s gisho P(Aw) d(Aw) ds 
+0 +0 
+. | | el(@—@o)s e7 isha P(Aw) d( Aw) ds. 


Because in the quasi-static approximation the variation in Aw is slow 
compared with the oscillations in the integrand, the contribution to the 
integrals arises solely from the region where the exponent is zero. 
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(Aw) o¢ P(Aw) 


as given by the simple argument above. 

We now wish to evaluate the occurrence probability P(Aw) in terms 
of the interatomic potential. This type of calculation was first performed 
by Holtsmark (1919) for the special case of broadening by charged 
particles (where the interaction is either the linear or the quadratic 
Stark effect) and by Margenau (1932) for the general case of scalarly 
additive interactions. We assume scalar additivity throughout this 
article for neutral atom interaction in absence of any knowledge of 
their non-additivity (apart from calculations with three interacting 
atoms by Dalgarno and Victor (1967)). 

Let us imagine a total of N; perturbers in a spherical vessel of 
radius Ro, centred about the radiating atom. The probability of find- 
ing a perturber in the spherical shell defined by R and R + dR is 47,R? 
dR3N7/(47 R3) = 3N,R? dR/R3. Thus we have that 


P(Aw) o HI ( 2 i R? ar,] (14) 
°K 


where R; is the position of ith perturber, and the integration is over 
these positions of configuration space which yield the correct pertur- 
bation Aw. The integral can, however, be written in such a form that it 
goes over all space with the use of a Dirichlet factor 


+ 0 
{= (x d(Aw)) Pips 
x 


—o 


Nr 


where Q) = : > V(R;) — Aw. This factor takes the value 7 when 


i=1 


Q = 0, i.e. when the total frequency displacement is Aw, and is zero 
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elsewhere; scalar additivity has been assumed at this point. We may 
therefore write 


Ro Ro 
Nr 

P(Aw) d(Aw) = (=) | af | R? R2... Ry, dR, dR, ...dRy, 

° 0 10) 

+00 ait 

| sin(x d(Aw)) wavgy 

x 
+0 Ro 


3 Nr Nr 
= (S| d(Aw) {{ [x dR errs} Ce aa 
. —o 


10) 


It is convenient to introduce the function 


Ro 
U(x) = fa — e!V(Rx/h) R2 dR (15) 
10) 


in which case we have 


+00 


P(Aw) = i (1 — 3U(x)/R3)%* e~ #42 dx, (16) 


Next we allow Nz and R, to tend to infinity, with the condition that 
N7/37Ro° = N. This is readily achieved by means of the expression 


Lt (1 — 3U(x)/R3)8? = Lt 0 — 420) NIN 7) = 6A, 


Nr>o Nr>o 
We therefore have 


+00 


P(Aw) = ieee e7 ixdo dx. 


—-7o 
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The integral may be expressed in a more convenient form by writing 
47U(x)N = A(x) + iB(x). (17) 
Then, recalling that U(— x) = U*(x), we obtain 


foe) 


T(Aw) oc P(Aw) = | e—4@ cos {B(x) + xAw} dx. (18) 


Oo 


Further progress depends on introducing an explicit form for the inter- 
atomic potential V(R), and this will be discussed in section 2C. 


B. Quantum Theory of Collision Broadening 


1. General formulation 

We turn now to an exposition of the quantum theory of collision 
broadening. It is important at the outset to understand the essential 
differences between the classical and quantum treatments of the prob- 
lem. We consider quantum effects in the perturbing and radiating atoms 
separately. The interaction between them is already in effect treated 
quantum mechanically in the classical theory, for V(R) can be the result 
of a quantum calculation, provided both states of the radiator can be 
considered stationary on the time scale of a collision. 

The first step is to examine whether the path of a perturber relative 
to a radiator can be adequately described by a classical orbit. The 
relevant parameter here is the ratio of the de Broglie wavelength of the 
atom to the minimum interatomic distance. This latter quantity is 
defined simply by the impact parameter for the critical ““Weisskopf”’ 
phase change. Closer collisions, of course, occur, but since they cause 
larger than critical phase changes in the radiation, the exact phase 
change and the details of the collision have no effect on the line profile. 
The condition for the classical path approximation to be valid for most 
collisions is, then, 


h h 
pee GmkTyt < Pe 
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(As we have already pointed out, this is also the condition which ensures 
that the path is straight.) Taking the worst case, we have to deal with 
experiments at a temperature of, say, 20 K with helium as the perturber. 
For this case the condition becomes 


pe > 3A. 


It is quite possible for this condition not to be very well satisfied, 
although most experiments are carried out at temperatures of several 
hundred Kelvin, in which case no difficulty arises. For complete genera- 
lity, however, the quantum description of the collision orbits must be 
used. 

Secondly, the question arises whether the use of the classical electron 
oscillator rather than a real atom to describe the emission or absorption 
of radiation is satisfactory. It is well known that this simple classical 
atom corresponds closely to a real atom with only two energy levels, 
or where the two energy levels involved in the transition concerned are 
completely uncoupled from all other levels. Except where there is 
degeneracy, or, what amounts to the same thing, where the collision 
itself causes transitions between energy levels of the radiator, we expect 
this approximation to be very good. Two kinds of degeneracy are 
concerned: one is the universal degeneracy associated with the quantum 
number m, and the other the much rarer accidental degeneracy where 
levels with different values of other quantum numbers have the same, 
or nearly the same, energy. If collision induced transitions occur, we 
expect that not only will the damping constant of the line be changed, 
but also the shape will no longer be Lorentzian. This is clear from the 
classical analogue, in which the collisions cause coupling between the 
oscillation concerned and other modes of oscillation of the system. 

The quantum mechanical analogue of eqn. (2) is 


Io) => lo = ay) CF 1d 1A>P iS (19) 


if 


where the initial and final states |i> and |f> are the time-dependent 
state functions referring to the whole system of radiators and perturbers. 
5(@ — w;,) is the Dirac 8-function, which takes the value 1 when 
w = wj;, and is zero elsewhere; p; is the number density of radiators in 
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the state |i>, and d the electric dipole moment. As in the case of classical 
theory, it is easier to calculate the autocorrelation function Qs). 
Formally, we can write, 


+00 


0s) = | #0) et joe Sites (<fldli> P pe (20) 
if 


It is convenient to introduce the time development operator, defined by 
Te) Gas (21) 


where H is the time-dependent Hamiltonian operator for the whole 
system. When T(t) acts on an eigenstate |a> of H we have 


IG) [a> = eo |'a> 


where £, is the eigenvalue of H associated with the state |a>. The 
autocorrelation function may therefore be written in the form 


96) = > <flaT@ |1><iH1dTO IP) Pi (22) 
if 


where the asterisk on an operator means Hermitian conjugate. 

No further progress can be made without introducing |i> and |f> 
explicitly. At this point we abandon a wholly general treatment, and 
invoke the classical path approximation. We refer to the article by 
Baranger (1962) for an account of the complete quantum mechanical 
theory: to our knowledge it has never been used in the interpretation 
of measurements of collision broadening by neutral perturbers. 


2. The classical path approximation 
We begin by writing the Hamiltonian for the complete system as 


H=H,+Hp+V 


where H,, Hp and V are respectively the Hamiltonians for the radiating 
atom, the perturbers and the interaction between the radiator and 
perturbers, and the wavefunction for the system as 


Pt) = x(t) $(0) (23) 
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where y(t) and ¢(t) are respectively the time-dependent wave functions 
for the radiator and the perturbers; y(t) depends on the coordinates of 
the radiator only and ¢(t) on those of the perturbers only. This wave 
function ignores the interaction between radiator and perturbers, and 
this influence must be evaluated by perturbation theory. The wave 
function describing the motion of the perturbers obeys the time depen- 
dent Schrédinger equation: 


ih as =H, 4(¢). 


To find the equation obeyed by y(t), we multiply eqn. (23) by 4*(t) and 
integrate over the perturber coordinates y,: 


= | HO 8 Osx 


from which we obtain 


h a = -{ W(t) H, $*(t) dx, + fou + Hy + V) ¥(t) $*() dx, 
= [H, + v(t)]}x(t) (24) 

where 
WO = | HOV 8) dx, (25) 


We can therefore proceed as though the atomic wave function y(t) is 
subject to a time varying perturbation given by v(t), and in the classical 
path approximation (that is with the d(t) assumed to be localized wave 
packets) v(t) is given simply by the classical potential describing the 
rapidly varying interaction between the radiator and perturbers. 

The effect of the perturbers may then be taken into account by writing 
a time development operator T, which operates on the perturbed 
atomic wave function y(t) and is defined by 


ih a 


= [H, + v(t)] T.@. 
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The autocorrelation function ®(s) now becomes 
P(s) = » [Ka |dTyts)|a><¢a|d Tits) |a> +c.0.], (26) 


where T,(t,s) = T,(t + s) — T,(t), a labels the upper states and a the 
lower states of the atomic transition, the average is taken over all 
times t and c.c. stands for complex conjugate. To clarify the analogy 
with the classical expression, we separate the time dependent from the 
time independent factors, and write ®(s) in the expanded form: 


®)= > Kald|a><B|d|b><b| Tas) |a> <a THO) |B> 
wd +c.c.], (27) 


If there is only a single upper state |a> and a single lower state |a> 
involved there is only one term in the summation and (27) reduces to 


P(s) = |<a|dl|a> |? [Ka| Tits) |a><a| TA8) |a> +0.c.], (28) 


which, apart from a numerical factor, is equivalent to eqn. (3). 
It is convenient for further development to work with the time 
evolution operator in an interaction representation. We write 


U(0,s) = exp(iH,s/h) T(0,s). 


This is precisely analogous to the device of working with Y(s)e~'®°* in 
the classical treatment. We then have 


ih a = exp ((H,s/h)v(t) exp(— iH4s/h)U = v'(t)U 


This equation may be solved by the standard perturbation expansion: 


8s s 


U(O,s) = 1 — (iff) | ¥(s1) ds, -+ (8)? | v'(s;) ds, { v'(s2) ds, +. 
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while in terms of U(s) the correlation function becomes 


06) = > Kaldla>< Bld] b> mean 


abap 


¢a | U*(4s) | B >< b/U(4s) | a> + cc.) (29) 


3. The impact approximation 

The calculation of the average in the impact approximation now 
follows very closely the classical calculation of section 2A 2. The 
change A@®(s) in ®(s) is written down; it involves A[U(t,s)], which is 
evaluated with the use of the impact approximation, namely that the 
two factors in A[U(¢,s)], = [U(¢,s)(U(t,s + As) — 1)], are statistically 
independent, and that the phase changes are caused by successive 
collisions. Just as in the classical case the line profile is Lorentzian 
provided there is no overlapping of lines, and the damping constant and 
shift are given by 


(oe) 


a oaND | (psce@ Comte needy (30) 
10) 
B = 2mNo | e-*-* sin (my — 15) pp G1) 
oO 
where 
ernst CaS a> (32a) 
et-im — <a | S| a> (32b) 
and 
+ 00 + 00 +0 
S = 1 — (i/h) | v'(s,) ds, + (i/A)? { v'(s1) ds, | v'(s2) ds. +... 


(33) 
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is the perturbation expansion of the S-matrix, defined by 
S = U(o, — oo) 


In these expressions, represents inelastic and 7 elastic (phase changing) 
collisions. If the inelastic collisions are ignored, ,. = 0, and the expres- 
sions for y and 8 reduce to those obtained with the use of classical 
theory. In this case the phase changes can be calculated without resort 
to perturbation theory. 

The neglect of inelastic collisions is certainly justifiable for foreign 
gas broadening where the energy levels of the radiating atom are so far 
separated that the frequencies present in the collision, of the order of 
the reciprocal of the duration of the collision, are smaller than those 
which correspond to the energy-level separations. Schuller and 
Oksengorn (1969b) have investigated the case of the alkali atoms with 
their fine structure doublets, but have concluded that the inelastic 
collisions cause a significant contribution to the broadening and shift 
only in the case of lithium. The fundamental reason for this is that, 
although the cross-section for inelastic collisions is quite large for 
transitions between certain fine structure components, it is, except for 
lithium, significantly smaller than that for phase changes of order unity 
(i.e. the Weisskopf cross-section). 

The case of self-broadening is quite different. Here the essence of the 
process is the exchange of the energy of excitation between the radiator 
and the perturber, and inelastic effects are all-important. A number of 
quantum mechanical calculations have been performed, the results of 
which are quoted in section 3B 2. 

The question of the effect of degeneracy with respect to the quantum 
number m requires some further comment. It arises even in the classical 
theory, where it takes the form of enquiring whether the radiating 
oscillator rotates through an angle z as the perturber flies past. If this 
always occurred, the result of every collision would be a phase change 
in the radiation of at least 7. Experimental evidence clearly rules this 
out, and it has usually been assumed that no rotation occurs for colli- 
sions which would otherwise produce phase changes of less than about 
unity. In quantum theory the question is raised by asking whether, in the 
laboratory coordinate system, the m-value of the radiator changes as a 
result of a collision. According to Baranger (1962), if the atomic 
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Hamiltonian H, is invariant under rotations and if the collisions occur 
isotropically, the question can be evaded, and the averaging over m 
carried out in a way which, in fact, amounts to what is usually done in 
classical theory: the interaction potential is simply averaged over m. 
The difficulty is, however, that this seems to assume that, in the co- 
ordinate system defined by the collision axis, a large number of tran- 
sitions between the several m-states takes place, and the question must 
be raised as to whether this is really so. The whole matter seems to be 
rather obscure, and further clarification is needed. 

If the collision broadened lines of a spectrum overlap, eqns. (30) and 
(31) for y and f are no longer valid. Off-diagonal elements of the 
S-matrix need to be calculated, and the matrix for the whole system of 
overlapping levels diagonalized. The line profiles are not Lorentzian, 
nor in general are they simply the sum of the Lorentzian profiles of each 
component in the overlapping system. Little detailed experimental or 
theoretical investigation of such cases has been carried out for neutral 
atom broadening; much more has been done for broadening by charged 
particles, but this falls outside the scope of this review. 


C. Interatomic Potential 


1. Introduction 

The theory of collision broadening has been described in sections 2A 
and 2B without reference to any specific interaction. We must now 
direct our attention to the choice of a suitable model for the interatomic 
potential. In the following paragraphs we make initially some general 
remarks about the interaction between pairs of atoms, both like and 
unlike, and then go on to describe the predictions of collision-broaden- 
ing theory in its various approximations when different potentials are 
assumed. , 

It was London (1930) who first evaluated by a quantum mechanical 
method, the long-range, attractive, van der Waals part of the potential 
between neutral atoms: the well-known dispersion contribution. He 
wrote the classical electrostatic interaction between two neutral atoms 
in a series expansion in R~* (where R is the internuclear distance) and 
treated this as a perturbation on the system. Since the results of this 
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calculation are rather different when the atoms are like and unlike, we 
shall treat the two cases separately. 

When the two atoms are unlike (as in foreign-gas broadening) 
London’s calculation shows that the first-order perturbation energy is 
zero, but that there is a second-order term 


apm Ret S faa fu 
2m? rT (Ey: — E,) (Ev — E)) (Ey: a ee yee eee 79 


lF-l’ 


(34) 


where one atom is in state |k>, the other in state |/>, f,,- is the oscil- 
lator strength of transition k —k’ for the first atom, fj, applies simi- 
larly for the second atom, and E,, E, are the energies of the states of 
the two atoms. (There are also higher order dispersion terms in R78, 
R-1'°....to which we shall refer below.) 

This expression may be greatly simplified if the perturbing atom has 
a much greater separation of its energy states than the optically active 
atom (that is, if it has a lower polarizability); this situation is frequently 
the case in collision-broadening experiments, where the spectral lines 
are often those of alkali or alkaline-earth elements and the perturbers 
are often rare-gas atoms. The term (E£,, + E,, — E, — E)) in the 
denominator may then be set approximately equal to (Z,, — E,), where 
I refers to the low-polarizability atom, and the summation separated 
into two summations, each involving only one atom. This procedure 
permits the expression of the second-order potential energy in terms of 
the dipole polarizability a, of the perturber, and we have 


AE=— Reo > |CkIrlk’ >? (35) 
k’ #k 


(see Uns6ld, 1955). In the case of a one-electron spectrum, this expres- 
sion may be further simplified to 


AE = — R-% e? ag ry? (36) 
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where r,2 is the mean squared radius of the orbit of the optical electron 


in state |k>. A value for r,2 which agrees with more sophisticated 
calculations to within 5%, is found using the formula 


ry? = 4a2 n*? {5n*? + 1 — 311 + 1)} (37) 


where dy is the radius of the first Bohr orbit, n* is the effective principal 
quantum number of the state |k> and / is the orbital angular momen- 
tum quantum number. In foreign-gas broadening studies, then, a 
suitable model for the long range interaction is 


V(R) = — CsR-° (38) 


where in many instances Cg = e” a, (r2? — r7) and the subscripts u and 
1 refer respectively to the upper and lower states of the transition. 

At short range, when the wave-functions of the interacting atoms 
overlap appreciably, a force of repulsion must become the dominant 
term in the interaction. It is convenient algebraically to assume that the 
potential in this range varies as R~+?, a form which has been used in 
kinetic theory as well as in the study of collision effects. That is, we write 


V(R) == - C12 Rane 


where C,, = C,2 (upper) — C;> (lower). A semi-empirical method of 
finding the constant C,, was suggested by Hindmarsh (1963), who made 
the assumption that the repulsive potential between any pair of atoms 
will be the same when the radial charge density betweeu them has a 
particular value chosen to give the best fit to the experimental results; 
that is 


Ci. = GRo- 


where g is a constant, Rp = R, + R, and R,, R, are the radii of the 
atoms 1, 2 where the radial charge density has the appropriate value. 
Hindmarsh found that using both Coulomb and Hartree-Fock wave- 
functions, a radial charge density of 0.012 a.u. to define R, and R, and 
a value for q of (0.9 + 0.3) x 10-73 J gave best agreement with 
experimental results. The degree of correlation between experiment and 
this simple theory, shown in Fig. 1, is remarkably good. 
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log (Cj2x10"'°) 
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Fic. 1. The correlation in experimental results between C,2 and R!* where 
R is an internuclear separation defined in the text. R is in and Ci2 
in erg A??. 


With these considerations in mind, it is clear that the Lennard-Jones 
(12,6) potential 


V(R) — Cs | a ee = Ce RS (39) 


forms a useful model, in terms of which foreign-gas broadening data 
may be interpreted. Indeed, this potential has been widely used in 
recent years since it became clear that many observed features could not 
be interpreted in terms of the single parameter van der Waals potential. 
The one major difficulty with the Lennard-Jones model has been the 
discrepancy between theoretical Cs; values and those experimentally 
derived: the Jatter generally exceed the former by factors of about 3. 
The possibility that this discrepancy may be caused by neglect of the 
higher order dispersion terms in R~®, R~'°... has been investigated 
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by Hindmarsh, du Plessis and Farr (1970), who used the (12,8,6) 
potential 


VR) = Op RY Cg Re = eee (40) 


and found that the situation was considerably improved in many cases. 

Another possible source of difficulty is the assumption of straight 
trajectories for the colliding atoms. This point has been investigated by 
Schuller and Vodar (1960a, b), who used classical paths in their calcula- 
tion, assuming a potential of the form 


V(R) => CO ca O, 
V(R) = — C, R-6, R>«. 


The effects of these refinements, together with the simpler potentials 
already described, are discussed below. 

A knowledge of the interaction between like neutral atoms is required 
for resonance broadening studies, where “resonance broadening” refers 
to the broadening of an excited atomic energy level through interaction 
with like atoms in the ground state, the two levels being linked by a 
strong dipole transition. (See, for example, Kuhn and Lewis (1969).) 

London’s perturbation calculation shows that the first-order term 
is not zero when the interacting atoms are of the same species; in fact, 
in many cases it is possible to neglect all terms of higher order for the 
purposes of pressure broadening. This leading term has the value 


; yerhf 


47mw9 


A= R= 


where R is the interatomic distance, f and w are the oscillator strength 
and angular frequency respectively of the resonance line. The factor y 
takes the value —2 if the magnetic quantum number m of the excited 
atom (relative to the interatomic axis) is 0, and the value +1 ifm + 1. 
Thus sum rules show that positive and negative values of AE occur with 
equal probability. The root mean square value over all m in the J 
representation has been given by Foley (1946), with the result that the 
J dependence of AE is expressed in the following equation: 


/3e7hf ( 25 1 le: 


AE = R~ a 
87mw, \2J’ + 1 


COLLISION BROADENING OF SPECTRAL LINES 169 


A proper understanding of resonance broadening requires quantum 
mechanical calculations, the results of which are discussed in section 
3B 2. 


2. The impact approximation 

With the assumptions of the impact approximation, classical collision- 
broadening theory, and also quantum theory if collision induced transi- 
tions are ignored, shows that the spectral line will always be Lorentzian 
in shape; the half-intensity width and shift, however, are sensitive to 
the choice of potential. Recalling the results of the calculations of 
section 2A 2, we have 


2y = 2Noo, (41) 
and 
B = Noo; (42) 


where o, and o; are given by eqns. (5) and (6). 

How the values of 2y and f depend on V(R) is examined below. 

(i) VR) = + C,R-’. For this general inverse-power potential, 
integration of eqn. (10) by parts yields the result 


eed) C, 1-—p 


= ———___— — 43 
n(p) = +7 TOD Ww” (43) 
from which it follows that 
o, = qrP/p—-A ee ee = ee ip = = *) cos ee 
ho I(p/2) p-—l p-—1 
(44) 
and 
Beer, ‘2a Rat gt a lee r(2 = ) a 7 
a ho I(p/2) p-—1 p—1 
(45) 


. providing that p>3. The sign of a, is always positive; the sign of o; is 
the same as that of V(R). Clearly then, the ratio |8/y| is tan (7/p — 1). 


| 
| 
| 
| 
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The connection between Weisskopf’s simplified treatment, which gives 
a zero shift, and the more detailed results of the Lindholm-—Foley treat- 
ment described here, may be understood readily. Weisskopf assumes 
only strong collisions to be effective (that is, collisions for which 7 > 1 
radian, and p > po, where 7(p,.) = 1) and neglects entirely the effect 
of weak collisions. Now for a strong collision, cos 7(p) and sin 7(p) 
are rapidly oscillating functions which average approximately to zero. 
Thus, strong collisions give 2y ~ 27Nip,? and B ~ 0. The weak, distant 
collisions, on the other hand, cause small phase shifts, so that cos 7(p) 
~1 and sin 7(p) ~ 7(p); from this it follows that 2y ~ 0 and B ~ 
27Nd<yn>, where <7n> is the mean phase shift occurring between two 
strong collisions. Close collisions, then, cause most of the broadening 
of the line and the distant collisions are responsible for the shift. 

(ii) V(R) = — C.R~°. Substitution of p = 6 and the assumption of 
the negative sign in (i) leads immediately to the expected broadening 
and shift for a van der Waals potential, so that 


2y = 8.16( Coin) oF rN. 
and 
B = — 2.96 (C,/h)?/> 09/5°N 


and the ratio 2y/8 = — 2.8. It was in fact the failure of many observa- 
tions to produce a ratio in even approximate agreement with this which 
led some investigators to explore the possibilities of the Lennard-Jones 
potential. 

(iii) VCR) = C,.R~'? — CgR~®. The Lindholm-Foley theory in 
the impact approximation was developed for the Lennard-Jones (12, 6) 
potential by Behmenburg (1964) and independently by Hindmarsh, 
Petford and Smith (1967). Here we have 


_ 630 Cae Sle 3a Ce HG -11 5) 
(0) me GaP 7 Sea a cee 


where 


ee (2 ss) 
~ PONS he (46) 
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and 


el 


PAS 
_- ee pas Soy as} 6/5 PAS: 
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_ following the notation of Hindmarsh et al. Thus, 


32\ 2/5 
2y = 89 (=) Nob?!> (C6/h)?!> B(a) 
and 
37r\ 2/5 
p25 ‘A Nb3/5(C,¢/h)2/> S(a) 
where 
B(a) = | x sin? {4(ax~11 — x-5)} dx 
10) 
and 


S(a) =| x sin (ax711 — x-5) dx 


18) 
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(47) 


(48) 


(49) 


(50) 


(51) 


The functions B(a) and S(a) are tabulated by these authors for a wide 
range of values of a; we show them here in graphical form in Fig. 2, 
along with the function $(a)/B(a), which is the ratio B/y. In order to 
interpret experimental results in terms of a (12, 6) potential, this latter 
curve is used with the experimental ratio to determine a, and hence 
S(a) and B(a). Then, with the use of either of the measured values B/N 
or y/N, Cg and C,, are found by a straightforward calculation. Unfor- 
tunately, the oscillatory character of the functions S and B sometimes 
leads to ambiguous results, although it is often possible to deduce the 


appropriate pair of constants. 
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B(a) 


Least logiol«) 


Fic. 2. The functions S(«), B(a) and S(a)/2B(a) computed for a Lennard- 
Jones potential (Hindmarsh ef al., 1967). 


(iv) V(R) = Cy,R~'? — CgR7® — CeR~°. The effect of the disper- 
sion term —CgR~® has been investigated by Hindmarsh, du Plessis and 
Farr (1970). The phase shift function becomes 


sr 63mCla! oe yee Geeta. 
1?) = aaa dost ole MIGGne ia 
Gy ex ee 
where a and x are defined in eqns. (47) and (46), and 


Sm (8\7/5 
=i (2) cxinerycu' (52) 
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Hence 
32\ 2/5 
2y = 8a (=) Nob?!> (C6 /h)?/> Bla, €) (53) 
and 
30 2/5 
p= 25 (=) Nb3!> (C6 /h)?/> S(a, &) (54) 
where 
Bia, é) = [> sin? 4 (ax711 — éx-7 — x75) dx (55) 
i¢) 
and 
S(a, €) = [> sin (ax—** — &x—7 — x77) dy. (56) 


8) 


Since it is clearly impossible to deduce all three constants from the two 
measured quantities y and 8, Hindmarsh et al. use a theoretical method 
for evaluating the parameter C,/CZ/° for each case investigated, using 
as a basis the work of Davison (1968), and requiring only the knowledge 
of the quantum numbers n* and / for the two states involved in the 
transition. Then for a given é, the functions S and B can be computed 
and the constants Cs and C,, deduced exactly as in (iii). 

(v) Numerically calculated potentials. Attempts have recently been 
made to carry out numerical calculations of the interatomic potential 
with sufficient accuracy to use in line-broadening theories. The methods 
used for evaluating the interaction of atoms which combine chemically 
are not sufficiently accurate: it should be noticed that the depth of the 
potential well in many interesting cases of foreign gas broadening is 
only a few meV or less. Two methods of calculation which aim to give 
sufficient accuracy in the relevant region of interatomic distance have 
recently been proposed. The first (Roueff, 1970) treats the short-range 
part of interaction as essentially arising from the scattering of the 
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valence electron of the radiating atom by the perturbing atom. The van 
der Waals interaction is added in separately. The expression for the 
short-range potential V,(R) involves the polarizability a, of the perturber, 
the binding energy U of the valence electron of the radiator, a scattering 
length L and the probability ¢?(R)dR that the valence electron is 
found in the spherical shell of radius R and thickness dR, where R is 
the distance between radiating and perturbing atoms: 


1 TragU 4a,U? 


f as TagU 
VR) =n oc {|} — tat ity) ee } 


where all quantities are expressed in atomic units. 

The second set of numerical calculations is due to Baylis (1969). 
Here the noble gas atom is again treated as a polarizable dipole in order 
to find the long-range interaction, which behaves asymptotically as 
R-°®. The short-range interaction is calculated by means of a pseudo- 
potential technique. The spin-orbit interaction is included, and so far 
only alkali-noble gas systems have been treated. 

Once the potentials have been evaluated, it is a comparatively simple 
matter to use the Lindholm—Foley theory to calculate the broadening 
and shift. At the time of writing this has been done for the potentials 
calculated by Roueff; comparison with experiment will be dealt with in 
section 3B 1. 


3. The quasi-static approxiimation 

In this approximation, described in section 2A 3, the intensity of 
radiation at any frequency is proportional to the probability of occur- 
rence of that frequency. Including the effect of multiple interactions 
(with the assumption of scalar additivity of the frequency perturbations) 
the result of these calculations is given by eqn. (18). We now apply this 
result to some specific interatomic potentials. 

(i) V(R) = + C,R-?. Provided that p > 3, we have, using eqns. 
(15) and (17) 

A(x) + iB(x) = SNe i (2 - 1] (C/Ay ee ( =| | x22 
3 Z Dp Dp 

a result obtained by partial integration. Further progress requires the 
choice of a specific value for p. 
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(ii) V(R) = — C.R-°®. The integral J(Aw) now takes the value 0 if Aw 
is positive (since all perturbations are tolow frequencies) and the value 


Qe 
I(Aw) = > (Colf)'! N | Seo |-9/? exp[— 4aCN?/9h | Av |] 
(57) 


if Aw is negative. 
The intensity maximum occurs at 


Vdd tan Cre 
Aw max - 7 jh aN Ss 
showing that the quasi-static shift is proportional to the square of the 
perturber density for a van der Waals interaction. (This contrasts with 
the linear dependence of the impact shift.) On the high frequency side 
of the maximum the intensity drop is predominantly exponential; in 
the far low-frequency wing the variation is proportional to |Aw|~3/?. 
This latter result was obtained by Kuhn (1934) who drew attention 
to the fact that the far wing of a line, which is produced by close colli- 
sions, is largely formed as a result of binary encounters. If then the 
frequency perturbation Aw = w — wo is caused by a single perturber 
at R, the intensity in the range w, w + dw may be equated to the 
probability of finding a lone perturber in the region R, R + dR. 


Thus I(Aw) dw oc R? dR x | Aw |~GOt?/? dw 


for a potential of the form C,R~?. 

(iii) V(R) = C,,R~12 — CgR~°. The more difficult problem of 
carrying out the quasi-static theory for a Lennard-Jones potential was 
first solved by Bergeon, Robin and Vodar (1952) and Bergeon (1954), 
who restricted their calculations to the region of high perturber density 
because of computational difficulties. Recently this work has been 
extended by Farr and Hindmarsh (1968) to the low density region. 

As before, this work is done with the assumption that the individual 
frequency perturbations may be added scalarly. It is convenient to 


P.Q,.E.—D 
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describe the potential in terms of its radius and depth parameters, o and 
e respectively, where 


Ges (C43/Ce)?” 
i ama (4C@/C, 2h), 


and then to work in terms of the reduced variables 


t? = 4ex (58) 
n = Aw/4e (59) 
and § = 7No3/3. (60) 


If we define a function 
A'(x) -- iB'(x) = [A(x) + iB(QX)]874, 


it is possible to write eqn. (18) in the form 
I(A)wo = [ te>°4") cos {dB'(t) -- nt*} dt (61) 


0 


apart from a constant factor. 
The function A(t) + iB’(t) is given according to Bergeon ef al. by 

the expression 
A'(t) + iB(t) = e'"/® T(— ) F(— 3 


1 
45 
Sa rah ges I(4) BE (oe oe ioe t it”) ee 
where F(a, b, z) is the confluent hypergeometric function 


az, aa+ by 2 


Flach ee ee 
(a, 8, 2) BLY Rb ey at 


whereas Hindmarsh ef al. use the equivalent form 


A'(t) de GBC) ee a - tk PAK — }) el-*/8-kr/4) (62) 


k=0 


i 
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The integral / (Aw) converges only slowly, especially at low values of 
perturber density (low 8), and can only be evaluated after a lengthy 
numerical computation, during which it is necessary to use an asymp- 
totic approximation for A and B at large values of t. The profiles com- 
puted by Bergeon ef a/. are shown in Fig. 3, and those of Hindmarsh 


High 


Low frequency 


frequency 


Fic. 3. The quasi-static theory intensity distribution for a Lennard-Jones 

potential. 6 = 1.5, 0.8, 0.3 in the three curves respectively. The highest 

density corresponds to about 100 atm and is the widest of the profiles 
(Bergeon et al., 1952). 


et al. in Fig. 4. The most noteworthy feature of this calculation is that at 
low 6, satellite bands appear on the low frequency wing at frequencies 
corresponding to wo — «, wo — 2, ... . They are, of course, here 
predicted to be much sharper than observed satellite features, as no 
impact broadening is included in the calculation (see section 2C 4). 
However, it is clear that the appearance of these bands, at Jeast in some 
cases, is a purely statistical phenomenon, arising from the fact that the 
potential goes through a minimum at an energy corresponding to a 
frequency of wo — «. 


4. Classical unified theory 
The difficulties inherent in a unified theory of pressure broadening, 
which describes the shift, width and shape of a spectral line under all 
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logiol L(y) ) 


Fic. 4. The quasi-static theory intensity distribution for a Lennard-Jones 

potential. 5 = 0.3, 0.1, 0.06, 0.01, 0.0014 for (i), (ii), (iii), (iv), (v) respec- 

tively. The lowest value corresponds to a pressure of somewhat less than 
an atmosphere (Hindmarsh and Farr, 1969). 


conditions and which is not subject to the restrictions of either the 
impact or quasi-static approximations, are such that no general analytic 
expression for the line shape can be found. The reader is referred to the 
review of Ch’en and Takeo (1957) for a description of the various 
approaches which had been tried before that date, all of which ulti- 
mately assumed V(R) = — C,R~°. Some recent attempts at approxi- 
mating the unified theory for a Lennard-Jones potential are given below. 
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Behmenburg (1968) has repeated the work of Lindholm (1945), 
replacing the van der Waals potential by V(R) = C,,R71* — C.R~°. 
Lindholm, whose work was the first unified theory to be proposed, 
writes 


[efn@-s)} = e7 Gs) + iK(s) (63) 


in which case the intensity distribution may be written, from eqns. (3) 
and (4) 


fo a] 


I(w) = | e~&) cos {K(s) + (w — wo)s} ds. (64) 


.0) 


The time-averaged exponential is then evaluated with the following 
assumptions: 

1. The atomic radiation suffers frequency changes, and therefore phase 
changes, as a result of collisions between radiators and perturbers. 

2. Perturbers, which are distributed randomly in space, follow straight 
paths with constant velocity relative to the optically active atoms. 

3. The duration of a collision is proportional to the impact parameter. 
4. Frequency perturbations caused by collisions which overlap in time 
are additive. 


5. The phase change suffered during a collision, [Aw(t)dt, is taken to 
0 


be a linear function of time. This is equivalent to approximating the 
complicated time dependence of Aw(t) by a constant frequency change 
during the collision. (See Fig. 6.) 

The method of evaluating [e'"”], is described in some detail by 
Lindholm, who shows that 


y | Pues 
G(s) = Noo, ja — cos n,) (s — 7,) — i sin 9, + 2r,| 
7) 
Tp=0 


D 


~ 2 
-/- : Noo, ja — cos «,S) (7, — s) — — sin «,s + 25| 
Kp 


ase (65) 


180 W. R. HINDMARSH AND JUDITH M. FARR 


and 


Ss ) | 
K(s) = > Noo, ja — cos n,) — + (s —7,) sin ns | 
Kp 
Tp=9 


2 : 
se > Noo, la — COS KS) — + (7, — 5) sin «5 | 
Kp 
ie 4 (66) 


where for a collision of impact parameter p, 7, is the collision time, 
x, the phase change per unit time during the collision, 7, = «,7, is the 
total phase change per collision and o, = 27pdp is the cross-section 
for such a collision. 

Now a collision with parameter p will produce (for an interaction 
of the Lennard-Jones type) the phase change given in section 2C 2 (iii) 
and, by assumption 3, 


T, = 2xp/b (67) 


where x is a constant of proportionality (~1) which Behmenburg 
determines later by adjusting theory to experiment, so that 


7 Cas 3 <*) 68 
x pies eerie a 
f 2yh \256<p'? 8 p® oe 
Using the parameters 
Pou: a) | 63 =, 
CAO bese (69) 
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Behmenburg substitutes these various quantities into the expressions 
given above for G(s) and K(s), to obtain for the intensity distribution 
in the line 


lo a) 


Hn) = | eG) cos {E'S hy + £915 AL») dy (72) 
0 
where hh = 4nNy(32C,/86)?/5 
and k= (w — wo) (x/6) (34C6/86)". 
The functions 4(¢, vy) and #(Z, y) are integrals given by 


HE” = | (If, GOO = 2) fe G2) ho G2) + 22) ede 


+ | ee 2 Ve — 7) "22 (0 2. a (6, 2) 29) 2dz 


y 


. 


WE, y) = | (1 — filS, 2) holS, z)) + f2(S y) (vy — z)] zdz 


+ [ (0 = 2G, 2.) hol 2) + a, 206 — y)] 2d 


where f7(.. z= cos (2029? —2z~!4)1,.46 2) =—sin 20277 — 271)] 
g, (¢, z,¥) = cos [2Ey(z—* — z~??)), 23(¢, z, y) 
= gin [(2£(z-* — z~**)] 
Bie Jig((; z) = 217/00 = 2°) 


| Only for sufficiently small y are the functions 4(%, y), ¥(Z, y) in closed 
mathematica! form, in which case they may be written 


H(f, y) = 0.450. €1/4 y1/4 + 0.43761/ y7/6 
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and ; 
wl, y) = — (0.187 £1/4 y/4 + 0.117 2/8 y7/®). 


For sufficiently large y, the functions are linear in y; at intermediate 
values they have to be computed numerically. 

The validity of this calculation may be tested by deduction of the 
special case of the impact and quasi-static limits and the far wing shape 
at low perturber density, and by comparison with experiment. The com- 
parison with experiment is described in section 3 of this review, and 
in the impact limit (small h and k) it may easily be shown that the 
intensity distribution here agrees with that of impact theory. To the 
authors’ knowledge the deduction from the general expression above 
of the quasi-static case and of the far wing shape has not been carried 
out. 

A different approach to the problem of the unified theory has been 
used by Fox and Jacobson (1969). Starting from the usual expression 


+0 


I(w) = | Peper ds 


for the intensity distribution, these authors write the correlation function 
(assuming scalar additivity of the perturbations): 


s 


®(s) i jexp {i i Aw(t) dt — 4 Mp mr} 


Avge 


where k is the Boltzmann constant, v? is the mean square relative 
velocity of perturbers and radiators and Me is the reduced mass. 
The averaging is then carried out over all initial positions and velocities 
of the perturbers, a procedure which makes the calculation essentially 
quasi-static in character. 

Fox and Jacobson introduce the time-dependence of the interaction 
by expanding the potential C,R~" — C,R~” in a power series in time 
(an idea suggested first by Holstein (1950)) and allowing the perturbers 
to follow classical paths with varying velocity during collision. Trun- 
cating the time expansion at terms in f°, they obtain for the correlation 
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function a complicated expression which has to be evaluated numeric- 
ally. A further numerical integration gives the line profile. An important 
feature of this calculation is the inclusion of the Boltzmann factor: this 
is important for close collisions, that is large interaction energies, and 
therefore for the extreme wings of the lines which are most easily 
observed at high perturber gas densities. 

This method of calculation produces the intensity distribution one 
would expect if all frequencies emitted had the same finite lifetime. This 
is clearly incorrect, as the configurations which give rise to these various 
frequencies have very different durations. At high perturber density this 
may well be of no consequence, since even a purely quasi-static calcu- 
lation predicts a smooth profile which will be relatively unaffected by 
the introduction of non-infinite lifetimes for the configurations. 

Takeo (1970) has made a calculation similar to that of Fox and 
Jacobson specifically for a Lennard-Jones potential, although the 
perturbers are assumed to follow straight paths. His calculated profiles, 
appropriate to perturber densities in excess of 10?°cm~%, show a 
satellite band on the low frequency wing, which however is in poor 
agreement with experimental profiles. 

A third different approach to the problem of finding a unified theory 
has been made by Hindmarsh and Farr (1969). Here the starting point 
is the quasi-static line profile, calculated for V(R) = C,,R7'? 
—C,R~°, which has already been discussed in section 2C 3. Using the 
fact that a Fourier analysis of a finite wavetrain of frequency w, 
averaged over random lifetimes, yields a Lorentz profile with damping 
constant y where y~* is the mean lifetime Hindmarsh and Farr give as 
an approximate expression for the line shape 


Fics yo’) /2 d dw’ 73 
i —e Cay erry ey bat 


where 7(w’)~! is the mean duration of wavetrains of frequency w’, 
and J,(w’) is the quasi-static intensity distribution. (The continuous 

change in w’ as a perturber flies past an optically active atom is not, of 
course, strictly equivalent to a sequence of values of w’ lasting for finite 
| durations.) 


P.Q.E.—E 


| 
| 
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For the purpose of calculating J(w) the line is divided into two main 
regions, the impact region where the quasi-static calculation is invalid, 
and the quasi-static region where /,(w’) forms a satisfactory starting 
point for the calculation of the profile. The boundary between these 
two regions on the low-frequency wing occurs at w,, where w, may be 
estimated as described in section 2A 2, eqn. (11). The contribution 
to I(w) from w > w, (note that w, < wo) is calculated by taking 
fI,(w’)dw’ within the impact region and giving it a Lorentzian distri- 
bution of width 2y and shift 6 calculated according to Hindmarsh 
et al. (1967) for the particular interaction potential under consideration. 
Since the quasi-static intensity on the high-frequency wing is very low, 
it may be set equal to zero, so that the line shape /(w) is given by the 
expression: 


(co!) yt Payee 
I(w u if w’ dw’ —s Sg oF een 
ahs eee Tomer Oo be tS |e eee 


(74) 


The mean lifetimes associated with the frequencies w’ < w, are esti- 
mated on the basis of an idea advanced by Kuhn and London (1934), 
namely that the time-varying frequency w(t) of acontinuously perturbed 
wavetrain may be considered to remain (0) for so long as the phase 
of w(t) does not differ from (0) by more than some arbitrary amount. 
Choosing this amount to be 1 radian, Hindmarsh and Farr define the 
semi-lifetime T of the frequency w(0) by the condition 


T 


[ eo de a0) = ie 


10) 


Then, replacing w(t) by the first non-zero term (term n) in the Taylor 
expansion of w(t) about the point t = 0, the condition for T becomes 


={ hea io. 
~ [(d"w(2)/d2"),=0 ; 


Using this expression, the mean lifetimes for all frequencies between 
w, and wo — e, where « is the depth of the Lennard-Jones potential, may 
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be estimated by considering only the attractive part of the potential; 
the mean lifetime for the frequency wy — « is calculated using the full 
expression for the potential, and this same value is used for all frequen- 
cies less than wy — « (further out on the low-frequency wing). This last 
approximation is not serious as the low value of J,(w’) in this region 


ensures that the major contribution to J(w) comes from frequencies 
greater than wo — e. 


-_ 


logiol 


-—!00 =(0) ! 


Wavenumber, cm 


Fic. 5. The quasi-static theory intensity distribution corrected for the finite 
velocity of perturbers. v = 2 x 10+, 5 x 10* and 10° cm sec~* for (i), (ii) 
and (iii) respectively. « = 60cm’, o = 6A (Hindmarsh and Farr, 1969). 


The final computation of J(w) is achieved by numerical integration 
of the expression given above. The comparison of this calculated line 
profile with an experimental profile is described in section 3; here 
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Fic. 6. The time dependence of the phase shift during a collision for a 
Lennard-Jones potential. The full line represents the true variation and the 
dotted line the linear approximation. 


(in Fig. 5) is shown the predicted behaviour of the low-frequency wing 
for a particular potential for three different values of 0. Clearly, the 
satellite band becomes narrower and more prominent at lower 
temperatures. 


3. EXPERIMENT 


A. General Considerations 


A great deal of experimental work in the field of neutral atom collision 
broadening has been undertaken in the 14 years since the appearance 
of the 1957 review of Ch’en and Takeo. The discussion in this section 
is not by any means intended to include an exhaustive list of these 
experiments. Rather, we have chosen a limited number which we con- 
sider to illustrate the way in which experiment has influenced the 
development of theory in this interval, the degree of accord between 
theory and experiment now attainable, the extent to which quantitative 
information on interaction potentials may be extracted from collision 
broadening measurements, and those phenomena which so far are 
unexplained. 


B. Impact Region Measurements 


1. Foreign gas broadening 

It is, of course, in the impact region where the experimental study of 
collision-broadened spectral lines is most easily interpreted quantita- 
tively. With the theories of section 2, it is possible to determine the 
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constants C, and C,, for an atom pair, where Co 12 = (Co6,12)excitea 
—(C6,12)srouna- The entire analysis rests on the assumption that the 
van der Waals, (12,6) or (12,8,6) potential chosen is a good model for 
the interaction between a perturber and an optically active atom in 
either state: having accepted this, a determination of the excited state 
constants is made possible by combining collision broadening data 
with that, for example, obtained from elastic scattering experiments 
which provide ground state values of the constants. 

It is usual for elastic scattering experiments to be interpreted in 
terms of a (12,6) potential, and a critical survey of available experi- 
mental results has been made by Bernstein and Muckerman (1967), 
who have published a list of ‘“‘best’”” parameters of the ground state 
potential for a number of atom pairs. Table 1 has been drawn up from 
these values, and shows in addition the calculated values of (C6) srouna 
found with Uns6ld’s formula, which is given in section 2C. 


TABLE 1. 
(CG) ground 
10m) 
€ Cc 
Atom pair 105729) 10-* cm exp theor 
Li-A 0.92 S57 20 18 
Li-Kr 1.5 3.6 26 28 
Li-Xe 2.4 3.8 45 44 
Na-A 1.0* 3.07 18* 20 
Na-Kr 1.4 4.5 80* 30 
Na-Xe 2.1 4.4 120* 49 
K-A 1.0 3.8 30 DT, 
K-Kr 1.4 4.5 90* 42 
K-Xe 1.8 <4 74 68 


* Less reliable values. 


Until about six years ago, in the absence of a theory of pressure 
broadening which coped with a more complex interaction, impact 
region experiments were interpreted exclusively in terms of the van der 
Waals potential, although as early as 1945 Kleman and Lindholm, 
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after an investigation of the behaviour of the sodium resonance lines 
under the influence of argon, had concluded that this interaction is 
inadequate at small distances. 

This conclusion was confirmed by Hindmarsh (1959, 1960) and 
Behmenburg and Kohn (1964). The latter authors measured the broad- 
ening of the resonance lines of sodium (37S,-3*P3,2) and strontium 
(51S -5!P,) by various gases under flame conditions at about 2500 K, 
and found, particularly with helium as perturbing gas, that the measured 
collision widths deviated significantly from theoretical values. 
Hindmarsh, also using helium as foreign gas, examined both the reson- 
ance and intercombination lines of calcium (41S,-4'P,, 44S -4°P1), 
which were both confined well within impact theory limits, and found 
that the ratio 2y/B, predicted to have the value —2.75 for a van der 
Waals interaction, was 35 + 15 for the resonance line and 7.5 + 2 
for the intercombination line. Concluding that short range repulsive 
forces must necessarily be taken into account for a meaningful inter- 
pretation of collision broadening, Behmenburg (1964) and Hindmarsh 
et al. (1967) published theoretical papers based on the (12,6) potential. 
These have been described in section 2C. 

Since 1964, several experimental studies of impact region collision 
broadening have been interpreted in terms of a (12,6) interaction. The 
first of these was Behmenburg’s (1964) study of the sodium line (37.$;2- 
37P3/2) perturbed by helium and by argon. In the case of argon, the 
shift-to-broadening ratio lies on the oscillating part of the S(a)/2B(a) 
curve, and so the interaction constants cannot be determined unam- 
biguously. In the case of helium, the measured C, of (33 + 10) x 
10~°° J cm® is anomalously large compared to the theoretical value of 
2.9 x 10-°° Jcm®, and the shift to a lower frequency is unusual for 
this perturber. However, the collision broadening of alkali-metal doub- 
lets is not well understood because of the possibility of collision-induced 
transitions, a problem which is discussed below. 

Singlet lines of alkaline-earth elements, being free of such effects, 
are more convincingly interpreted, and Smith (1967) investigated the 
calcium intercombination line (41S)-4°P,) formed in absorption in 
atmospheres of helium, neon, argon and krypton at temperatures of 
approximately 1030K and perturber densities up to 7 x 1018 cm73. 
In all cases, impact region conditions were well satisfied with shift and 
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broadening varying linearly with perturber density. The results were 
interpreted using the Hindmarsh et al. (1967) version of the Lindholm- 
Foley theory for a (12,6) potential. The following interaction constants 
were obtained: 


6 
Atom pair LOmoce pom ounml Omaozel Cina 
Ca-He 2.0 aa 0.2 0.97 
Ca-Ne 1.8 + 0.2 Let 
Ca-A 10.8 + 0.8 TA 
Ca-Kr 25 ae 2 54 


These Cg values are all higher than the theoretical values by a factor 
of about 2, a result which is attributed by Hindmarsh et al. (1970) to 
the neglect of higher order dispersion terms in the interatomic inter- 
action (see below). 

A searching test of the validity of the (12,6) potential in the interpre- 
tation of collision broadening was applied by Vaughan and Smith 
(1968), who studied the red line of krypton (4p°(?P3/2) 5s[3]. — 4p° 
(?P3/2) Sp[$]2) perturbed by low pressures of helium, neon, argon 
and krypton at temperatures of 80K and 295K. They were thus able to 
test the complex temperature dependence predicted by Hindmarsh 
et al. as well as to measure the interaction constants. The shift-to- 
broadening ratios at the two temperatures were fitted as closely as 
possible to the theoretical S(a)/2B(a) curve, with the restriction that 
each pair should show the correct temperature dependence. The 
corresponding values of S(a) and B(a) were then used to predict the 
temperature effects in shift and broadening individually. These are 
compared in the following list: 


Shift at 295K/Shift at 80K Width at 295K/Width at 80K 


Perturber Observed Predicted Observed Predicted 
He 2.62 2.16 1.87 ies 
Ne 0.93 0.84 1.67 1.70 
A 2.29 1.79 2.07 1.64 


Kr 1.47 1.43 1.48 1.43 
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Only in the case of argon does a serious discrepancy arise, which 
Vaughan and Smith suggest may be due either to the inadequacy of the 
theoretical treatment for potentials which are neither deeply attractive 
nor sharply repulsive, or to the wrong choice of interaction potential. 
To test the latter possibility, they calculated the S(a), B(a) curves for 
(10,6) and (18,6) potentials, but neither made much difference in the 
argon region. The observed C, values are in satisfactory agreement with 
the Unsdld calculations: 


Ce, (10- °° J cm®) 


Atom pair Observed Predicted 
Kr-He 4.5 4.7 
Kr-Ne 15.9 8.6 
Kr-A 24.5 35.4 
Kr-Kr 54.6 53.8 


The study of alkali-metal doublets has been pursued since the early 
days of experimental collision broadening (Fiichtbauer and Gédssler, 
1933; Watson and Margenau, 1933; Hull, 1936; Peterman, 1933; 
Ny and Ch’en, 1937), and a striking feature of all such measurements 
is the existence of apparently unpredictable differences between com- 
ponents of the same doublet perturbed by the same gas. Schuller and 
Oksengorn (1969a) have calculated the j-dependence of the broadening 
and shift in the impact approximation (using only dispersion terms in the 
interaction) and have shown that the 7P,,. component behaves as in 
Lindholm’s theory, but for the ?P3,. component there is a variation due 
to the anisotropy of the excited state interaction. Applying their 
calculation to the case of the broadening of the sodium resonance lines 
by argon, they find that the calculated shifts are in good agreement with 
those measured by Kleman and Lindholm, but that there is a discre- 
pancy in the half-widths. The interference with the collision broadening 
process by collisional excitation transfer between the fine structure 
levels is regarded by some authors as being fundamental to this problem, 
and both theoretical (Callaway and Bauer, 1965) and experimental 
(Krause, 1966; Jordan and Franken, 1966) evidence shows that the 
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cross-sections for these two processes are of comparable size. However, 
Schuller and Oksengorn (1969b) show that this process has a marginal 
effect, measurable only in lithium. 

A recent comprehensive study of the pressure broadening of the first 
few lines of the principal series of caesium by helium, neon, argon, 
krypton, xenon and carbon tetrafluoride has been undertaken by 
Ch’en, Garrett, Gilbert and co-workers (1966 I, Il; 1967 III, IV; 1969 
V, VI). Measurements have been made over a wide density range, from 
the impact region up to pressures of 300 atm in some cases. The general 
features of these results are explained by recent classical unified theories 
(see below). The low density measurements, the results of which are 
summarized in Table 2, along with the results of some of the earlier 
work mentioned above, are not interpreted quantitatively because of 
some unusual features, which include non-linearity of the broadening vs. 
density curves near the origin and an apparently very strong temperature 
effect on the broadening at higher densities. There is a possibility that 
both of these effects are experimental: the non-linearity arising from 
the uncorrected effect of hyperfine splitting and instrumental slit width, 
which amount to some 0.6cm~', and the temperature effect arising 
from analysis of lines whose peak absorption is in excess of 50%. There 
is some evidence that heavily absorbed profiles can cause overestimated 
line widths (McCartan, 1970). 

Rostas and Lemaire (1971) have measured the caesium doublet 
6*S—7?P broadened by helium and argon. These measurements, made 
at low pressure, represent an improvement on those of Ch’en and 
Garrett, and are quoted in Table S. 

McCartan and Hindmarsh (to be published) have made a detailed 
investigation of the second doublet of the potassium principal series 
(4781/2 — 5?P1)2, 3/2) formed in absorption in atmospheres of helium, 
argon and krypton under strict impact region conditions. The results 


TABLE 3. 
Helium Argon Krypton 
Width Shift Width Shift Width Shift 
4044 P3;4 7.08 + 0.2 | 0.37 + 0.02] 8.0 +0.1 |—2.22 + 0.03] 5.50 + 0.15 ]|—1.97 + 0.05 
4047 Bs 4.94 + 0.04 | 1.22 + 0.02 | 7.54 + 0.05 |—1.75 + 0.08] 6.6 + 0.2 1.84 + 0.05 


Units of 10-2° cm-! (cm-3)-1, 
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of this study are summarized in Table 3. It was found that not all of 
these could be fitted to the S(a)/2B(a) curves, a fact which the authors 
attribute to the inadequacy of the “straight path’? approximation in 
these cases. The interaction constants deduced by McCartan et al. are 


Ce J cm® 
Pi)2 P32 
He “= Owes (SC= 
A — 6.6% 107° 
Kr 4.3 x 10=*% AS Ono 


Ce J cm ~ 
Prj2 P3)2 
— 1:6 A105 20° 
—_ Sex Ree 
12S 100° DEX MOmTeo 


Combining these figures with the ground state measurements quoted 


earlier, the excited state constants are 


(GAR Iie 
P32 P3)2 
He = 9.5 x 10-5 
A — 6.9 x 10-% 
Kr 52x 10-% 5.5 x 10-6 


(Cia)ex J Gin? 


Pij2 P3)2 

— WO OAL 

—— Sie Se 1K Be 
4 S<10822 OSea1Omiwe 


The experimental facts regarding collision broadening in alkali 
doublets may be summarized as follows: 


1. The difference in shift of the fine structure components depends 
on series member and perturber, not on the nature of the alkali. 
2. Where there is a high-frequency shift, the *P,,. component is 
shifted more than the 7P3,/,. The ratio is larger for the first series 


member than the second. 


3. Differences in shifts to low frequencies are usually small, and there 
is apparently no way of knowing which component shows the 


stronger shift. 
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4. The broadening of the 7P3,. component is normally slightly 
greater than that of the ?P,,. component, although exceptions 
are found. 


An important consideration in assessing the usefulness of collision 
broadening measurements is the degree of accord between the measured 
constants and the theoretically expected values. The constant C;, is 
not a reliable guide, as the “theoretical” value is in fact a semi-empirical 
value (see section 2C). The constant Cg is a different matter, as theor- 
etical estimates are more reliable, and any discrepancy must throw doubt 
on the validity of collision broadening theory. Table 4 shows, in 
columns four and five, the sort of agreement which is found between 
theoretical C; values and those deduced from collision broadening 
experiments using the Lindholm—Foley treatment for a (12,6) potential. 
On the whole, the experimental C, values exceed the theoretical ones 
by factors of about two or three. The last column shows the experi- 
mental values when a (12,8,6) potential is assumed. There is in general 
an improvement, which lends support to the idea that higher terms in 


TABLE 4. 


Ces experimental 


10mcoicmes 
C; theoretical 

Author Line Perturber | 10-°° J cm® (12,6) (12,8,6) 
(1) Ca 6573 He 0.062 0.2 + 0.02 0.10 + 0.02 
(1) Ca 6573 Ne 0.113 0.18 + 0.02 0.10 + 0.02 
(1) Ca 6573 A 0.46 1.08 + 0.08 0.23 + 0.02 
(1) Ca 6573 Kr 0.70 2.53 + 0.2 0.61 + 0.06 
(2) Ca 4227 He 0.30 0.3 +0.1 0.18 + 0.06 
(3) A 6965 A 2:4 9 +4 2.9 +3 
(4) K 4047| Kr 38.1 Sy? +15 30 +9 
(5) Sr 4607 A 2.62 43 +04 2.8 +0.8 
(6) Kr 7601 He 0.47 0.45 + 0.09 0.36 + 0.07 
(6) Kr 7601 Ne 0.86 1.59 + 0.3 0.81 + 0.15 
(6) Kr 7601 A 3.54 2.45 + 0.5 1.05 + 0.2 
(6) Kr 7601 Kr 5.38 5.46 + 1.1 44 +1 


(1) Smith (1967); (2) Hindmarsh (1959); 
(4) McCartan and Hindmarsh (1969); 
(6) Vaughan and Smith (1968). 


(3) Hindmarsh and Thomas (1961); 
(5) Farr and Hindmarsh (1971); 
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the dispersion interaction (R78, R~'°...) may not be negligible in the 
excited state. Finally, the calculations of Roueff (1970) have been com- 
pared with experiment by Rostas and Lemaire (1971) for the caesium 
line A4593 (671). — 77P1/2) broadened by helium and argon. The 
calculations, which do not include spin-orbit interaction, are not 
expected to apply to the other component of the doublet, 44555, which 
was, however, also measured. Table 5 compares the results of Rostas 
and Lemaire with the calculation of Roueff, and shows the agreement 
to be highly satisfactory. 


2. Resonance broadening 


Resonance broadening calculations in the impact approximation 
predict that pure resonance broadened lines are Lorentzian in shape 
with half-intensity width in radian s~! of 


oFa=1 ie 4nef,,,.N 


2y = kj; (7 Aer 


MW o 


where N is the atom density, f;,,and w, the oscillator strength and angu- 
lar frequency respectively of the resonance line and /,j’ are the angular 
momentum quantum numbers of the ground state and resonance state 
respectively. The width is independent of temperature, and the shift 
is predicted to be zero. (See the review of resonance broadening by 
Kuhn and Lewis (1969).) The value of the constant k,;, depends on 
j,j and the approximations made in the calculation. Some recent 
values are: 


ki)2, 3/2 

ko,1 k1)2, 3/2 kij2, 1/2 Ki 2, 112 
| Aliand Griem (1965, 1966) 1 1 1 1.00 
| Omont and Meunier (1968) 1.04 = = pad 
1.04 1.04 0.903 1S 


| Stacey and Cooper (1969) 


' Because of the difficulties of self-absorption in emission lines, most 
| measurements made directly on resonance lines (largely in the case of 


| 
| 
| 
| 
| 
| 
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TABLE 5. 
Mo on ee ee Ee ee 
Experiment Theory 
Rostas and Lemaire Roueff 
Cs line Gas (1971) (1970) 
4593 2y/N 8.8 + 0.5 9.2 
67S4)2-77P 1/2 He B/N +1.5 + 0.1 
2y/B +5.85 +6.3 
A4555 2y/N 6.9 + 0.4 
67S )2-77P3/2 He B/N +0.73 = 0.05 
2y/B +9.45 
A4593 2v/N 6.7 + 0.5 4.9 
67S1j;2-77Pij2 | A BIN —1.63 + 0.05 
2y/B —4.1 —11.6 
A4555 2v/N 5.8 = 0.3 
67.S1)2-77P3/2 | A B/N —1.55 + 0.06 
2y/B —3.74 


y/N and B/N are in units of 10-7° cm~+/cm~?. 


alkali metals) have been made in absorption. Moser and Schulz (1959) 
have tabulated the experimental results available at that time, and 
discrepancies amounting to several orders of magnitude occur between 
different experiments and between experiment and theory. Since then, 
Ya’akobi (1966, 1968) has made measurements on the resonance lines 
of lithium, but these results also seem to be not very quantitative. 

A more accurate technique involves the indirect study of the broaden- 
ing of a resonance level by making measurements on an emission line 
whose lower state is the resonance level. Quantitative results were first 
obtained this way by Hindmarsh and Thomas (1961), who investigated 
the self-broadening in argon and were able to estimate the oscillator 
strength of the resonance line A1048 by studying the self-broadening 
of the line 48264. Lewis, Rebbeck and Vaughan (1969) adopted the 
same method and measured the resonance broadening of the 4?2P 
states of potassium by studying in emission the lines (47P,/2,3/2 
—n? S12, = 6, 7, 8). This gave the values k4)2,3,..= 1.33 + 0.3, 
Kyj2,1;2 = 1.24 + 0.3, with the ratio k1)2,3/2/k1)2,1;2 = 1.09 + 0.06. 
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Kuhn and his co-workers (Kuhn and Vaughan, 1964; Kuhn and 
Lewis, 1967; Vaughan, 1966a, b, 1968; Stacey and Vaughan, 1964) have 
made a major investigation of resonance broadening in the inert gases 
helium, neon, argon and krypton using this indirect technique, with 

the object of testing the results of theory. In all cases studied the varia- 
tion of half-width with density was very close to linear over a density 
range of 10'© — 4 x 1017 cm~. The results obtained for lines leading 
to the upper state of a strong resonance line are given in Table 6. 


TABLE 6. 
Element ws Transition 2y/N at 80K 2y/N at 280K BIN 
A 10-29 cm-1(cm™ 3)—4| 10-29 cm-1(cm-3)- 1 | 10-20 cm-1(cm-3)-1 

He 7281 | 21P,-31S, 7.85 + 0.4 8.1 + 0.6 0+ 0.1 

6678 | 21P,-31D, 7.8 +0.4 Thess USS) 0 + 0.05 
Ne 6678 |1P,-?P, 6.27 + 0.15 6.2 + 0.3 OF Oe 
A 7503 1P,—*P, 15.66 = 0.7 15.0 + 0.7 OF O02 

7514 | 3P,-?P; 3.6 +0.3 3.9 + 0.4 0.07 + 0.03 
Kr 7685 | 1P,—-5p, 10.8 + 0.5 11.0 + 0.6 0+ 0.5 

7587 | 3P,—5p 12.8 +0.2 13.1 + 0.6 0.1 + 0.1 


In all cases there is large broadening, no measurable shift and no 
measurable temperature effect. Other lines were measured which are 
transitions to a weak resonance level, and these show 5-10 times smaller 
broadening, an increase of broadening with temperature and a non- 
zero shift: in these lines a higher order interaction interferes with the 
resonance broadening. 

The case of helium is particularly interesting as the oscillator strength 
of the resonance line has been computed very accurately by Schiff and 
Pekeris (1964), so that the measured resonance broadening leads to a 

value of ky,,;. This value, 0.99 + 0.06, is in excellent agreement with 
recent theoretical values; Omont and Meunier (1968) considers the 
value 1.45 to be accurate to + 0.04. Using this value of ko,, oscillator 
strengths of the resonance lines of neon, argon and krypton have 
_ been derived from resonance broadening measurements, and are in satis- 
factory agreement with theory. 
An interesting phenomenon has shown up in these rare gas experi- 
_ments, namely that the Lorentz width extrapolated to zero density is 
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significantly higher than the expected natural width. This “extrapola- 
tion anomaly” has been explained by Kuhn, Lewis and Vaughan (1965) 
in terms of an additional coupling width. 

The technique of resonance fluorescence may be used to make accurate 
measurements of self-broadening in certain favourable cases. Such ex- 
periments involving *P, resonance levels have been performed by Omont 
and Meunier (1968) and by Happer and Saloman (1967). The k,;, values 
are different for the orientation and alignment signals; the theoretical 
value for the alignment broadening is k = 1.32. Omont’s study of the 
6°P, level in mercury, and Happer and Saloman’s investigation of the 
3P, resonance in lead both give good agreement with this value of k 
within experimental accuracy, and the latter authors further find the 
ratio of the broadening constants for the orientation and alignment 
signals to be 1.21 + 0.05, compared with a theoretical ratio of 1.25. 

Little work on resonance broadening has been done at higher densi- 
ties where the impact theory is invalid. The quantum mechanical 
calculations of Mead (1960) and Reck, Takabe and Mead (1965) predict 
a lowering of the resonance level and a linear variation of y with 
/N. These facts are not inconsistent with the experiments of Vaughan 
(1968) and Lauriston and Welsh (1951). 


C. Non-impact Region Measurements. 


Since Minkowski (1935) first showed that the intensity distribution 
in the far low-frequency wing of the sodium resonance lines perturbed 
by argon was of the Aw~ */? form, numerous experiments have provided 
data concerning the non-impact region of pressure broadening: we 
discuss here only a few of the more recent of these. 

Granier and his co-workers have presented, in a series of papers (1961, 
1963a, b, 1965a, 1966a, b), the results of a study of the influence of 
perturber mass and temperature on the broadening and shift of the 
resonance lines of mercury and the alkali metals at high perturber 
density, and have come to the conclusion that, especially when dealing 
with light atoms, repulsive forces cannot be neglected. This result, of 
course, parallels the work of Behmenburg and Hindmarsh in the impact 
region. Proper comparison of theory and experiment outside the impact 
region, however, has only been possible since the development of a 
(12,6) unified theory, and the experiments which have been quantita- 
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tively assessed in this way are so far only those of Ch’en et al. on the 
behaviour of the caesium resonance lines when perturbed by argon and 
neon (see section 3B), and the study of McCartan e¢ al. of the low 
frequency wing of the potassium line (47S )-5*P/2) perturbed by 
krypton. 

Ch’en’s work has been interpreted by Behmenburg (1968), who has 
considered the effects of both neon and argon, and by Fox and Jacobson 
(1969), whose calculations apply to broadening by argon only. Figures 
7, 8 and 9 show respectively the half-intensity width, shift and asym- 
metry of the caesium line (67S)-6*P,/.) as a function of relative density 
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Fic. 7. A comparison of theory and experiment in the pressure ee nee 
of Cs A8943 by argon (Behmenburg, 1968). 
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of argon (density relative to that at N.T.P.) obtained experimentally 
by Ch’en and calculated by Behmenburg. The C, and C,, values used 
in the calculation are obtained by Behmenburg from the low-density 
measurements (corrected for the effects of hyperfine structure and slit 
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Fic. 8. A comparison of theory and experiment in the pressure shift of Cs 
A8943 by argon (Behmenburg, 1968). 


width) using the well-established (12,6) impact theory. The value of the 
factor x (see section 2C), which fixes the assumed collision time, is 
chosen by Behmenburg by forcing a fit between the experimental and 
theoretical parameters at some characteristic point. In the case of argon, 
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Fic. 9. The asymmetry of Cs A8943 as a function of the argon pressure. 
The inner scale refers to the van der Waals potential (Behmenburg, 1968). 


this point is that at which the asymmetry is a maximum. The calcula- 
tions are thus based on the values 


C= 5.0 <A0s** Dem? 
C5 — 3.0 x 10> 48° J cm® 
vy == 1.01: 


Fox and Jacobson obtain, for the same line, the theoretical curves 
shown in Figs. 10 and 11 (along with the experimental points) by using 


Ce = 2.25 x 10-5 Jom 
Cro = 4.59 x 107199 J cm?? 


The theoretical value of C, is 3.6 x 10~°° Jcm°®. 

Behmenburg has also interpreted the neon broadening, fixing y by 
forcing a fit between experiment and theory where the relationship 
between shift and density goes through zero for the second time. The 
results, shown in Figs. 12, 13 and 14, are less good than for argon, 
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Fic. 10. A comparison of theory and experiment in the pressure broadening 
of Cs A8943 by argon (Fox and Jacobson, 1969). 
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Fic. 11. A comparison of theory and experiment in the pressure shift of 
Cs A8943 by argon (Fox and Jacobson, 1969). 
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Fic. 12. A comparison of theory and experiment in the pressure broadening 
of Cs A8943 by neon (Behmenburg, 1968). 


which Behmenburg feels may be due to the inadequacy of the (12,6) 
potential in describing the neon-caesium interaction. The constants 
used are 


Cra 59x. 10-"" Jem? 
C12 = 942 x 10749? J cm? 
x = 1.18 


and the theoretical value of C, is 0.91 « 10~°° Jcm®, 
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Fic. 13. A comparison of theory and experiment in the pressure shift of 
Cs 48943 by neon. The inadequacy of the van der Waals potential is 
particularly evident in this case (Behmenburg, 1968). 


McCartan et al. have, in contrast to Ch’en, studied a low-frequency 
wing profile at low perturber density, but at such large displacements 
from the unperturbed line that the impact theory is no longer applicable. 
Using the quasi-unified theory of Hindmarsh et al., they have inter- 
preted the shape of the (47S,/.-5*P,/2) absorption line of potassium 
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Fic. 14. The asymmetry of Cs 8943 as a function of the neon pressure. The 
upper scale corresponds to the van der Waals potential (Behmenburg, 1968). 


broadened by krypton as arising from a (12,6) interaction potential 
where 


a2 fe lem ~ 
o = (8.0 + 0.3) x 10-8 cm, 
which is in excellent agreement with the impact region result that 
e=25+2cm7! 
= (0.0 = 0,9). 107* cm, 


Figure 15 shows that the calculated and experimental profiles are in 
very good agreement over an intensity range of at least four orders of 
magnitude, including the prominent satellite band at about 25 cm? 
from the core. The contribution to the profile beyond the band is 
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almost entirely due to “impact” broadening of the band itself, due to the 
finite duration of the wavetrains absorbed at that frequency. There is 
some indication of a feature at about 50 cm~*, which would be expected 
on the basis of the quasi-static calculation of Fart et al. (1968) which 
predicts satellites at 2e, 3«, .. 

This work has removed any serious doubt as to the nature of the 
alkali-rare gas low-frequency bands, of which a great number of 
observations have been made and several theories advanced. The 
essential experimental facts about such bands, all of which are consistent 
with the theory of Hindmarsh et al., are as follows: 


1. Satellites are more likely to be observed the higher the polariz- 
ability of the perturbers (i.e. the larger the value of a,). 

2. At high pressures the bands become diffuse and disappear. 

3. The position of the satellite relative to the main line is almost 
independent of pressure and temperature (see, for example, 
Granier and Granier, 1966c). 

4. The bands are usually several wavenumbers wide even at low 
pressures, but the intensity drops sharply on the low-frequency 
side of the satellite. In general, the more distant the band the more 
diffuse it is. 


Recent studies of the low-frequency bands in alkalis include those of 
Ch’en and Wilson (1961), Ch’en and Fountain (1964), Takeo and 
Ch’en (1964), Jefimenko and Wiliams (1965) and Gwinn, Thomas and 
Kielkopf (1968), all of which provide data which are consistent with the 
idea of associating a band with a turning point in the potential well. The 
low-frequency bands observed by Ch’en, Atwood and Warnock (1961) 
on lines of indium and thallium perturbed by xenon and krypton also 
appear to be bands of this type. It is not, however, clear that the low- 
frequency satellite bands which have been observed near forbidden lines 
of alkalis in the presence of xenon are so easily explained. The first 
such observation was made by Lapp (1966), who reports an intense 
band in association with the weak caesium absorption line (6S-SD). 
This same band has been observed in emission by Gwinn, Thomas and 
Kielkopf (1968), and a similar band near the rubidium absorption line 
(5S-4D) (Besombes, Granier and Granier, 1969). No adequate theor- 
etical treatment of this phenomenon is available. 
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The remaining feature of foreign gas-broadened spectral lines which 
is to be discussed is the appearance, usually at high perturber density, 
of satellite bands on the high-frequency wing of a collision-broadened 
line. Most work on this topic has been carried out on the resonance 
lines of xenon, mercury and the alkali metals. The main features of 
these bands are: 


1. They are fainter than low-frequency bands. 

2. The heavier the perturber, the closer the high-frequency band is 
to the main line, although in general they are much further from 
the line centre than low-frequency bands. 

3. For a given atom pair, the high-frequency band separation is less 
for higher members of the principal series, and low series members 
usually show a second, fainter band. 

4. Structure is occasionally observed (Robin, 1957). 

5. The position of the high-frequency band appears to be temperature 
and pressure dependent (Besombes et al., 1969). 


Ch’en and Wilson (1961) and Ch’en, Atwood and Warnock (1961) 
give a comprehensive survey of band positions. 

Michels et al. (1956, 1957, 1959a, b) have studied the shape of the 
absorption line of mercury at A2537 (1.S)-3P,) under the influence of 
helium, neon, argon and krypton at pressures of several hundred 
atmospheres. In every case, features are found in the high-frequency 
wing. With helium as perturbing gas, a shoulder develops at 225 cm7? 
from the line when the relative density exceeds 150. This is attributed 
by the authors to a locally increased absorption at frequencies corre- 
sponding to the minimum in the ground state (12,6) potential, where 
the radial distribution of perturbers has a maximum. Neon has a similar 
effect on the line, producing at pressures in excess of 200 atm a pro- 
nounced high-frequency shoulder at 235 cm7}. 

With argon and krypton, the situation is rather different. Argon 
produces two high-frequency satellites at 86 and 170cm~?, krypton 
three at 34, 73 and 128 cm7?. In both cases, the whole pattern shifts 
to a lower frequency in the same way as the pressure is increased, 
maintaining approximately constant peak-to-peak separations. These 
results have been interpreted by Michels, de Kluiver and Ten Seldam 
(1959b) in terms of weakly bound mercury-inert gas molecules. They 
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estimate the contribution to the absorption coefficient from such mole- 
cules to be considerable in some cases, and find that this model is 
compatible with the measured variation of satellite intensity with 
temperature and perturber density. Taking a suitable potential, the 
high-frequency wings of the outer satellites of mercury/argon and 
mercury/krypton are described as a result of transitions from the lowest 
vibrational level of the ground state to the repulsive part of the upper 
state. 

Granier and Granier (1965b) have also studied the mercury 2537A/ 
helium high-frequency satellite, and find that the satellite separation 
(at low densities approximately 125 cm~*) increases with temperature 
at constant density, and with density at constant temperature. These 
results are not inconsistent with the ideas of Leycuras (1966) that the 
high-frequency satellites are formed as the result of adding to the line 
frequency a frequency Avp, which is the oscillation frequency due to 
the local perturbation of the charge cloud of the perturber by the optical 
electron of the active atom. 

The resonance line of xenon, 41469, perturbed by hydrogen and rare 
gases has been investigated by Granier, Castex, Granier and Romand 
(1967) in absorption, and in emission at low density by Rupin and 
Robin (1964). Their results give the satellite band positions in cm™' as 


Granier et al. Rupin and Robin 
Xe-H, 170; 350 — 
Xe-He 160; 300 not observed (density too low) 
Xe-Ne 140; 270 162° 255 
Xe-A 105; 250 102; 246 
Xe-Kr — 42; 171 


A tentative explanation of these bands has been offered by Marteau, 
Granier, Vu and Vodar (1967), who suggest that they may arise by 
combination transitions between the electronic resonance line and 
translational bands. In support of this they offer measurements of the 
frequency maxima of translational bands for the pairs Xe—He 
(160 cm~*), Xe-Ne (125cm~*) and Xe-A (65cm~*), although no 
firm conclusion is possible as these latter results were obtained at 
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pressures at least one hundred times greater than the pressures used in 
the collision broadening experiments. 

Most of the work on high-frequency satellite bands has been with the 
alkalimetals. Robin (1957), for example, studied the effect of high pres- 
sures of argon and nitrogen on the resonance lines of potassium, rubi- 
dium and caesium, and found that there is very little difference between 
the effect of the two perturbers, although nitrogen, the lighter gas, pro- 
duces rather more distant satellites. In each case the satellite relative 
to the line appeared to move initially to a low frequency, then back 
again as the pressure changed. 

Ch’en and Wilson (1961) have attempted an explanation of the 
systematics of the high-frequency bands in the principal series of 
lithium perturbed by helium, using the idea of an adiabatically split 
repulsive potential as proposed by Klein and Margenau (1959). The 
positions which they estimate seem to be satisfactory, but the intensity 
is incorrect. Jefimenko and Williams (1965) found high-frequency 
satellites associated with each of the first three potassium doublets 
formed in absorption in atmospheres of foreign gas. The resonance 
line bands were very diffuse; those at the second and third doublets 
much more intense. Jefimenko et al. attribute these bands to peculiari- 
ties in the potential curves during collision, an explanation which 
Granier and Granier (1966d) also feel is likely to be correct for the 
rubidium—neon high-frequency bands which they investigated. Their 
evidence suggests that each of the two components of the rubidium 
resonance doublet has associated with it one of the two observed high- 
frequency satellites. Most principal series doublets, however, seem to 
give only a single high-frequency satellite, although structure is some- 
times observed in the satellite. On the other hand, Ch’en, Smith and 
Takeo (1960) report that in sharp series each component of the doublet 
has its own high-frequency satellite. 

High-frequency satellites associated with the first forbidden doublet 
of rubidium (5S-4D) at 5167A have been observed by Besombes, 
Granier and Granier (1970). Three high-frequency satellites are observed 
with helium as perturber and two with each of the other rare gases. The 
high intensity of these satellites compared to the line itself leads 
Besombes et al. to feel that, as with “forbidden”’ low-frequency satellites, 
a simple explanation in terms of potential curves is not sufficient. 
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The theoretical understanding of the high-frequency satellite bands 
is so far very incomplete. The possibility that, like the low-frequency 
bands, they arise from turning points in the difference potential cannot 
be ruled out. The difficulty is that rather little is known about the 
detailed shape of the interatomic potential at the close distances which 
must be responsible for the formation of high-frequency bands. 

It is appropriate that this review should end on a note of uncertainty. 
Much has been discovered over the past three-quarters of a century; 
but the phenomenon of pressure broadening continues to bring to the 
surface problems of real significance for our understanding of the inter- 
actions of atoms with each other and with electromagnetic radiation. 
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